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A model of sequential resonant tunneling transport between two-dimensional subbands that takes 
into account explicitly elastic scattering is investigated. It is compared to transport measurements 
performed on quantum cascade lasers where resonant tunneling processes are known to be domi- 
nating. Excellent agreement is found between experiment and theory over a large range of current, 
temperature and device structures. 



Resonant tunneling in semiconductor heterostructures has motivated a large number of experimental and theoretical 
studies. One of the most studied case is the resonant tunneling diode, where a quantum well is formed by a double- 
barrier regioEpl Under a variable applied electric field, the transparency of the system is probed by coherent tunneling 
of electrons at the Fermi energy of the contact regiorPEl. The voltage-current curve exhibits a clear maximum when 
the emitter electrode aligns with a resonance of the well. This perfect case might be realised when the exit barrier 
is made so thin that the escape tunneling rate is much faster than other dephasing mechanisms. Because of strong 
in-plane scattering this condition is difficult to achieve. Usually current proceed by sequential tunneling as it is the 
case in quantum cascade lasers. 

In the pioneer work of Kazarinov, the current is expressed in a density matrix model where the resonance curve is 
found lorentzian with an homogeneous broadening given by the average value of elastic scattering matrix elements. 
As a consequence of the averaging the electrons tunnel between subbands conserving their in-plane wavevectors. 

More recently a refined model that include previously averaged-out second-order mechanisms was developecP. 
Second order scattering is known to yield gain without a net population inversionPI through scattering assisted optical 
transitions, but it also affects more generally resonant tunneling, by allowing transitions between subband states of 
different wavevector^. It is found that resonant tunneling occurs with conservation of the energy rather than the 
wavevector, contrarily to the first order casd^. 

In this letter we demonstrate the important role played by second order mechanisms in sequential resonant tunneling 
and therefore in the carrier transport of semiconductor heterostructures. 

When second order terms are considered in the calculation, the current density between a pair of subbands coupled 
through a barrier is expressed aS 
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with q± = h^2rn* (e^ ± A), where f\, is the carrier distribution in subband i at wavevector fc, A is the detuning 
between the subband edges A = €2 — ei, hSl is the coupling energy through the barrier, d is the difference between 
the two centroids of the wavefunctions d = Z2 — zi, e is the elementary charge and 7^ is the broadening of state i at 
wavevector fc. 

When a low density of electrons is distributed thermally in each subband, with the same electronic temperature T 
and that a same and uniform scattering potential is considered (7^ = 7fc =7)1 the current density can be integrated 
and simply rewritten as: 



+0(-A)(e-'3''l^ln2-ni)} 

where d{x) is the Heaviside function, with 0{x~) ~ 0, 9{x^) = 1 and 6{0) = \, (i — 1/kT with k the Boltzmann 
constant and Ui is the net population of subband i. 

The current density is no more driven by the population difference 712 — ni but by an effective population term. We 
want to examine two extreme case: equally populated subbands n2 = ni = n and one empty subband ni — 0. The 
first case is shown in Figjlja). The current density is dispersive shaped around the resonance. A negative current 
peak occurs when the detuning is negative, this is when the edge of subband 1 is above the edge of subband 2. When 
the subbands are aligned the current is zero and the first order approximation is recovered. The current then turns to 
be positive, after the edge of subband 2 has overcome the edge of subband 1, this is when the detuning is positive. The 
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dispersive shape is the consequence of electron tunneling at a constant energy rather than at a constant wavevector. 
As shown in FigQa) the first order model yield a zero current for any detuning. This case illustrate a superlattice: 
the current is zero until second-order scattering terms have been taken into account. 

The case where one subband is empty is shown in Fig|ljb). For negative detunings the current between the 
subbands is exponentially reduced as only the electrons with a sufficient kinetic energy are able to tunnel to subband 
1. Contrarily, for positive detunings, the first and second order curves overlap perfectly as all electrons are above the 
edge of subband 1 . 

Generally the first order model is recovered as the thermal energy largely overcomes the detuning energy kT ^ ^|A|, 
as it is the case irP. In Eqj2] the exponential cut-off tends to 1 as the temperature tends to infinity, spreading electrons 
in an uniform distribution. 
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FIG. 1: (a) Effects of second order contributions on tunneiing between a pair of equaiiy popuiated subbands. (A) Wiien tlie 
detuning is negative, the subband edge of subband 1 is above the edge of subband 2. As tunneiing conserves energy, the current 
flow from 1 to 2 is greater than the current flow from 2 to 1, yielding a negative net current. (B) When the subbands are 
aligned, the detuning is zero and both contributions cancels, yielding a zero net current. (C) The detuning is positive and 
therefore the edge of subband 2 is above the edge of subband 1, yielding a positive net current between subbands. (b) Empty 
subband 1. (A) The current is reduced as only a fraction of electrons can tunnel. (B) Models overlap perfectly. 



As second order mechanisms strongly affect the resonant current between a pair of subbands, we aim to show its 
impact on more complex semiconductor heterostructures like quantum cascade lasers-^' We therefore have implemented 
second-order effects in the computation of the voltage-current characteristic. 

The computational model is based on the density matrix where dissipation is included as rate equations for the 
populations and as dcphasing times for the polarisations. The precise implementation will be detailed somewhere 
els^. A typical quantum cascade laser is a repetition of a fundamental period as shown in Figj2] These periods 
are coupled through an injection barrier. Electrons are injected by sequential resonant tunneling from a period 
to the next one. The period itself can be separated in an active region where the laser transition occurs and an 
injector region where carriers are relaxed before they are injected into the next period. For many structures, and 
the ones presented here, the active region is coupled to the injector region through an extraction barrier as shown 
in Fig|2] We therefore have implemented tight-binding and sequential resonant tunneling at the injection and at 
the extraction barrier. In the active region an d the injection region, the carriers are relaxed through intcrsubband 
scattering. The mechanisms we have considered^^lli] include LO-Phonons, interface roughness and ionized impurities 
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(dopants) scattering. Non-parabolicity effects and self-consistency of the potential are accounted by the model. A 
uniform electronic temperature is computed for all subbands, based on the electron energy balanc^^^. The numerical 
simulations output the populations of the subbands and the current density flowing through the heterostructure. We 
are therefore able to predict the voltage-current characteristic of a particular quantum cascade structure. In order 
to test the impact of second order transport, we have implemented both first and second order sequential resonant 
tunneling models. 
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FIG. 2: 

Each structure is shown at injection resonance field. The layer sequence starts from the injection barrier and the thicknesses 

are in nm; roman, resp. bold, numbers indicate Ino.53Gao.47As, resp. Alo.48lno.52As alloy, acting as well, resp. barrier 
material, (a) Layers: 4.3/1.7/0.9/5.4/1.1/5.3/1.2/4.7/2.2/4.3/1.5/3.8/1.6/3.4/1.8/3.0/2.1/ 2.8/2.5/2. 7/3.2/2. 7/3.6/2. 5. 
Underlined layers are 1.5 x 10^^ cm~^ Si doped. Nominal sheet carrier density is 1.2 x 10^^ cm~^. Period length is 68.3 nm, 

repeated 35 times. The optical transition occurs at «154 meV. (b) Layers: 
4.8/3. 6/0. 2/3. 6/3. 5/5. 1/1. 1/5. 0/1. 2/4. 5/1. 3/3. 5/1. 5/M/ 1-6/3.3/1.8/ 3.2/2.1/3.0/2.5/3.0/2.9/2.9. Underlined layers are 
3 X 10 cm~ Si doped. Nominal sheet carrier density is 3.03 x lO" cm"^. Period length is 68.6 nm, repeated 35 times. The 

optical transition occurs at «167 meV. 

We present two quantum cascade structures in different coupling regimes. The first show in FigMa) has a strong 
coupling between the active and the injector regions as the extraction barrier is made sufficiently thin (22 A) . Contrarily 
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FIG. 3: 

The current curves are shown both in log scale (left axis) for inspection of low currents and in linear scale (right axis) for 
inspection of the dynamic range, (a) Experimental data (full line), simulation with second order resonant tunneling (dotted 
line), simulation with first order model (dashed line), (b) Single quantum well current-voltage curves for three temperatures: 
80 K, 180 K and 300 K. Measurements displayed in full line. Simulations in broken line. 



The current- voltage curves for both structures are shown in FigjS] The measurements are taken in continuous mode 
for low currents and in pulsed mode when current flow causes a heating of the sample. The cryostat temperature for 
the first sample is 300 K, while it is 80 K, 180 K and 300 K for the second. 

In Fig|3ja) the second order voltage-current curve fits the experimental data much better than the first order 
approximation from the very low currents to the maximal current. 

If we focus on the low field values of the voltage-current curve, the second order model clearly better predicts the 
experimental behaviour than the first order model does. In particular it yields a zero net current at zero field which 
is an important validation of the computational model. 

The result for the single quantum well structure are shown in Figj3f^b). We have shown simulated curves with 
second order model only, because the first order model failed to converge at low fields values and is largely off from 
the measurements. Apart from a constant serial resistance (0.5 fl) than higher systematically the experimental bias, 
the model was able to reproduce nicely the experiment, in particular for high temperature where the transport in the 
structure is clearly dominated by optical phonons. The model predicts the low temperature curve with less accuracy 
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because the transport at such a temperature also require the computation of scattering rates due to acoustical phonons 
and electron-electron interactions, which are not computed in the present model. 

Agreement between computed and experimental current- voltage characteristics has already been reported for other 
model approachs such as based on Monte Carlc*^ or non-equilibrium Green's function^^^ the comparison was done 
however on a much more limited range of currents, temperature and structure design. 

Formally, the current driven by tunneling between two subbands through a barrier or by optical absorption are 
physically equivalent because both processes conserve the in-plane wavevector. As a result, the striking agreement 
between the predictions of the second-order model and the experiment can be interpreted as a strong experimental 
evidence for the validity of the Bloch gain model. 
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